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Abstract: Air showers induced by cosmic rays create nanosecond pulses detectable at radio 
frequencies. These pulses have been measured successfully in the past few years at the LOw- 
Frequency ARray (LOFAR) and are used to study the properties of cosmic rays. For a complete 
understanding of this phenomenon and the underlying physical processes, an absolute calibration 
of the detecting antenna system is needed. We present three approaches that were used to check 
and improve the antenna model of LOFAR and to provide an absolute calibration of the whole 
system for air shower measurements. Two methods are based on calibrated reference sources and 
one on a calibration approach using the diffuse radio emission of the Galaxy, optimized for short 
data-sets. An accuracy of 19% in amplitude is reached. The absolute calibration is also compared 
to predictions from air shower simulations. These results are used to set an absolute energy scale 
for air shower measurements and can be used as a basis for an absolute scale for the measurement 
of astronomical transients with LOFAR. 

Keywords: Antennas; Large detector systems for particle and astroparticle physics; Spectral 
responses. 
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1. Introduction 


When a high-energy cosmic ray hits the top of the Earth’s atmosphere, a cascade of secondary 
particles is created that propagates through the atmosphere to the ground. This phenomenon is 
called an extensive air shower. A promising method for the detection of air showers is based on their 
radio emission. Air showers emit short broad-band pulses that are caused by the electromagnetic 
component of the showers [|I]]. As the shower develops in the magnetic field of the Earth, a time- 
varying transverse current of electrons and positrons is created, which causes radio emission [0. 
Also, air showers accumulate electrons, resulting in a charge excess in the shower front. This again 
induces a time-varying current, which leads to emission [Q]. The emission from both currents adds 
up coherently in frequencies below 200 MHz that correspond to the dimensions of the shower front, 
i.e. a couple of meters 0-0. Measuring the radio emission of air showers has shown to be a good 
method to determine the type and the energy of the primary particle e,®. 

The EOw-Erequency ARray (EOEAR) is used to study the radio emission processes of air 
showers le-imi- An understanding of the electric field of the radio pulse in each antenna is essen¬ 
tial for the reconstruction of the air showers. Eor this, the response of the electronic system and 
antenna to the incoming signal is most important. The recorded voltage traces depend severely on 
the instrumental chai'acteristics, such as the antenna gain, filter amplifiers or losses in cables. In 
addition, the gain-patterns of the used antennas are direction dependent, which has to be considered 
for the different arrival directions of the cosmic rays. In order to fully understand the phenomenon 
of radio emission and the underlying physical processes, an absolute calibration of the measured 
field strength is needed. This also enables us to test the absolute energy scale of different experi¬ 
ments such as EOPES Tunka-Rex [1T311 , or the Auger Engineering Radio Array (AERA) [jT^l, 
[1T311 of the Pierre Auger Observatory [|T0. Such tests are independent of air shower simulations 
and allow us to compare the results of the calibrations conducted by different groups. 

Eor astronomical observations of sources, essentially at infinite distances, a number of calibra¬ 
tion procedures are used by the EOEAR Collaboration. They mostly depend on known point-like 
astronomical sources which are resolved in imaging and used for flux calibration. Thereby, they 
are based on integrated or correlated quantities, requiring a significant amount of data. The signals 
of air showers are non-repeating pulses of a duration of less than 100 ns and usually just 2.1 /is 
of data are recorded per antenna and measurement. Therefore, a calibration based on resolved 
point-sources is challenging with our data. 

We present the results of two measurement campaigns that were conducted at EOEAR in order 
to understand and calibrate the low-band antenna together with the EOEAR system (Section @ and 
0. Both include reference sources that were placed in the vicinity of and measured with EOEAR. 
One method is used to verify the antenna model with respect to the directional as well as frequency 
dependent response of the antenna (Section ||). The other is used to provide an absolute calibration 
of the system. This is accompanied by a third method that is based on calculations using all-sky 
emission models of the Galaxy (Section 0. The results are also compared to those obtained with 
the same method at other experiments uni-iiM]- Eor reference purposes also a comparison to the 
absolute field strength as predicted by air shower simulations is shown in Section 0 
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Figure 1. Left; The six innermost stations of the Superterp. The black crosses indicate the positions of the 
LBA antennas, the gray full squares indicate the LORA scintillation detectors and the open squares mark 
the HBA tiles. One circular cluster of LBA antennas and two clusters of HBA tiles belong to every station. 
Right: A low-band antenna at LOFAR. The supporting vertical PVC pipe is shown with the LNA on top of 
it. The dipole wires are located to the right and left of the LBA and are attached to the wired ground-plate 


2. The LOFAR antennas 

LOFAR antenna stations are distributed over the Netherlands and Europe 0, Nine inter¬ 
national and thirty-eight Dutch stations make up an interferometric array of dipole antennas with 
effectively all-sky coverage. The core of LOFAR, located in the North of the Netherlands, con¬ 
sists of 24 antenna stations within an area of 12 km^. The Superterp consists of the six innermost 
stations inside the core within a diameter of 320 m as shown in Figure |I]. 

A LOFAR station comprises two types of antennas, the Low-Band Antennas (LBAs) [^TH , as 
well as the High-Band Antennas (HBAs) [P^H . Together, they cover the frequency range from 10 
to 240 MHz. Every core station in the Netherlands consists of 96 LBAs, 48 HBAs and 96 digital 
receiver units (RCUs) 0. The LBAs are grouped into an inner and outer circle that each contain 
48 dipoles and are called LBA Inner and LBA Outer, respectively. All read-out electronics are 
equipped with RAM buffers on transient buffer boards (TBBs) that store up to 5 s of data. These 
buffers can be read out after a trigger or on manual request. At the Superterp, an array of 20 
particle detectors is used to trigger a read-out of the LOFAR TBBs whenever an air shower above 
a minimum energy threshold has been detected [^31] . The data are then available as individual 
time-domain signals per antenna and are combined in one file for every station. 

2.1 The Low-Band Antennas 

The LBAs are designed to measure in the frequency range between 10 and 90 MHz 0, [^TH , where 
the lower end is determined by the strong ionospheric reflection of terrestrial radio-frequency inter¬ 
ference (RFI), and where the higher end is marked by the commercial FM radio band. Practically, 
the operational range is therefore limited to 30 — 80 MHz. 
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An LBA in the field of LOFAR is shown in Figure |T]. The sensitive elements of the antenna are 
four copper wires that act as slanted dipoles. The antenna is held upright with a PVC pipe and with 
tension on the springs in which the copper wires terminate. Two coaxial cables that carry the signal 
of the two polarizations run through the PVC pipe. These cables are also used for the power supply 
of the low-noise amplifier (LNA) that is on top of the LBA and acts as a voltage amplifier, including 
a notch filter for the FM band [^TJ]. According to the length of the wires, the resonance frequency 
of the LBA dipole should be 58 MHz. However, this does not include the influences of fhe ofher 
electronics. A different resonance frequency is obtained when including the whole analogue and 
digital system in the model. The antennas are known to be sky-noise dominated. Observations of 
CAS A and 3C295 were used to confirm that the fraction of sky-noise over total noise is above 0.5 
for most of the LBA frequency range [0,0. 

Every LBA has two polarizations. One polarization of the antenna is aligned from northeast 
to southwest, while the other is aligned orthogonally from northwest to southeast. The omnidi¬ 
rectional response, which is provided by the two orthogonal polarizations, allows for simultaneous 
observations of the whole visible sky. 

2.2 System characteristics and antenna simulation 

In order to obtain an absolute calibration for the LOFAR system, a frequency and direction depen¬ 
dent factor is needed that converts the measured voltage traces (in uncalibrated units) to the electric 
field emitted from the incoming air shower. Two contributions have to be considered: the vector ef¬ 
fective length of the antenna (VEL), and the additional frequency dependent gain and phase factors 
introduced by filters and cables.’ The VEE contains all directional sensitivity, as well as a large 
part of the frequency dependent gain and group delay. Understanding the antenna will therefore 
allow us to compare signals of air showers arriving from different directions. In order to achieve 
this, the entire system has to be considered. 

Every incoming electric field signal E{t) can be separafed into its two components perpendic¬ 
ular to the direction of propagation, indicating its two polarization components in eg and : 

E{t)=Ee{t)ee+E^{t)e^. (2.1) 

Eor the influence of the VEE of the antenna H that converts these electric fields E to voltages V, 
one can also consider them in the frequency domain rather than in the time-domain^ 

V{v,e,^)=H{v,d,^)-E{v,e,^), (2.2) 

where 6 denotes the zenith angle and <p the azimuth angle as indicated in Eigure @ and V is the 
frequency. 

Eor each frequency the VEE can be described by the Jones matrix J with its components for 
the X- and the T-dipoles of the EOEAR EBA 

^ (t: t:) ■ ft) ■ 

'This section summarizes briefly some aspects of antenna theory and introduces our notation. More detailed descrip¬ 
tions can be found in e.g. OSS 

^For this analysis the focus is on the amplitude response. Future work will also include the phase response of the 
system to transient signals. 
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Figure 2. Spherical reference coordinate system of a low-band antenna of LOFAR. Point T indicates the 
position of the transmitting source which is used for the gain calibration measurements of the LBA. The 
distance between transmitting and receiving antenna is denoted by r. Symbols 0 and (j) indicate the zenith 
angle and azimuth angle, respectively. 




Figure 3. Example response of the Vector Effective Length, in the form of output voltage AV for an LBA. 
Left; Response to an incoming wave polarized in the eg direction (circles), and a wave polarized in the ^ 
direction (triangles), with arrival direction 0 = 45° and (j) = 65°. Right; \Jxe \ component at 60 MHz as a 
function of direction for a wave fully polarized in the eg direction {Eg = 1). 


Since both dipoles are independent, the equation can be separated into the two components. For 
LOFAR, the Jones matrix of the LBA is modeled with the electromagnetic simulation software 
WIPL-D mTll and a customized software model of the gains of the electronics. This model contains 
some simplifications, but describes the system to our current best knowledge. 

The model predicts the response to an incoming wave with an electric field strength of 1 V/m, 
and is computed with steps of 1 MHz in frequency, 5 ° in 0 and 10 ° in (/). The components at inter¬ 
mediate values are obtained via trilinear interpolation when needed. An example of the directional 
sensitivity of the Jones matrix, as well as the frequency dependence are shown in Figure ^ 
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Figure 4. Equivalent circuit of the LBA. The circuit represents a voltage source including an internal resis¬ 
tance which is equal to the antenna impedance. The LNA, which works as an operational amplifier, is located 
directly behind the antenna. The equivalent circuit consists of a capacitor C and a resistor R. The output 
resistor Rout has to be matched to the impedance of the LNA, Zout- This matched situation corresponds to 
the calibration as performed with the flying reference source. 

We can visualize the antenna system of the LBA together with the load impedance of the LNA 
in an equivalent circuit which is illustrated in Figure ^ Vemf represents the incoming signal, the 
voltage at the inputs of the LNA is denoted by Vin, the antenna impedance by Zq and Gy is the 
voltage gain. As long as the input impedance Z,„ is large compared to the antenna impedance Z^, 
there will be no voltage drop across the antenna impedance. The ratio of the output and the input 
voltage can be calculated by 



(2.4) 


with 


(2.5) 


l+j(0RC' 


and ft) = 27rv. The load impedance of the LNA Zout and the output resistance Rom are 75 Q. for 
the LBA, respectively. This is also the typical resistor of the coaxial cables used in the LOFAR 
stations. Hence, the losses due to reflection are minimized. 

3. Calibration set-ups 

For the analysis of all measurements and to describe the radiation pattern of an antenna, a spherical 
coordinate system is used as shown in Figure ^ The origin of the coordinate system is located 
below the center of the antenna and all quantities refer to the electromagnetic far-field region. The 
emitted signal coming from large distances can be described as a plane wave. The vector of the 
electromagnetic field has an eg and etj, component in the coordinate system of the transmitting 
antenna. The signal propagates in the er direction. 
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Figure 5. Schematic view of the experimental set-up to measure the antenna gain pattern in a closed loop. 
The grey half circle denotes the path of the transmitting antenna. The angle 9 indicates the zenith angle of the 
transmitting antenna with respect to the LBA. The black dashed-dotted line indicates the cable connection 
between the LBA and the amplifier used for the calibration of the network analyser. 

It has been checked that all measurements have been conducted in the far-held region. The 
required minimum distance arises from far-held condition within the LBA operation bandwidth 
from 30 — 80 MHz and given the antenna dimensions D = 2.76 m to 

2D2 

r > —^ Ri 4 m. (3.1) 

A 

Consequently, with measurements taken at a distance of at least 12 m away from the source, the far- 
held approximations are valid for the measurements. However, the far-held criterion also includes 
the condition that A < < r to avoid spill-over from the near-held into the far-held regime. For this 
condition the distances might be on the small side. It has been shown before that the spill-over 
effects were not observed for this particular experimental set-up DTEH , while they were still relevant 
for other measurements HTTH . In this analysis they will be treated as additional uncertainties. 

In the two experimental set-ups, reference transmitters are used to calibrate the absolute scale 
of the LBAs. The emitted signals are signihcantly stronger than any contribution from system 
or ambient noise. Still, the signal strength was carefully adapted to not saturate the system and 
stay well below the non-linear regime of the amplihers. All measurements were performed using 
regular LBAs installed in the LOFAR core. 

3.1 Measurement of the antenna gain pattern 

The determination of the gain pattern of the LBA requires a study of the direction- and frequency- 
dependent sensitivity of the antenna to an incoming signal. Measurements between a transmitting 
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Figure 6. Drone and reference antenna as used at LOFAR. Left; Octocopter drone with the transmitting 
antenna mounted below. Right; Biconical antenna BBOC 9217 used for the calibration measurements. 

reference antenna attached to a flying drone and a receiving LBA at different zenith angles have 
been performed. More details about the analysis of these measurements can be found in Chapters 
4 and 5 of ra. 

3.1.1 Experimental set-up 

The in-situ calibration measurements of the LBA have been carried out at two antennas of station 
CS302 in the LOFAR core during a period of two days in September 2012. The set-up requires 
an LBA as receiving antenna and a transmitting antenna aligned in parallel, as well as a vector 
network analyzer. A schematic view of the set-up to measure the zenith angle dependence of the 
antenna response is depicted in Figure ^ 

For the measurement, a Rohde & Schwarz FSH4 vector network analyzer with two ports was 
used, where one port was attached to a bias tee that supplied the receiving antenna with 6 V DC. The 
other cable was directly connected to the LNA of the LBA and the received signal was measured 
with the vector network analyzer. Via the second port a signal was sent through the cable to 
the transmitting antenna attached to the drone. This signal was amplified (25 dB) in order to be 
significantly stronger than the background noise. 

In order to calibrate the network analyzer, the footpoints of both antennas were connected di¬ 
rectly. This was done only once before the actual measurements, and is represented by the dashed- 
dotted line in Figure ^ Thereby, losses due to the electronics of the calibration set-up and additional 
group delay caused by the cables are directly taken into account during measurements. The vector 
network analyzer has a systematic uncertainty on the gain of 0.6 dB. 

For the actual measurements, the transmitting antenna has been located above the LBA with 
the help of an octocopter drone, with the antenna mounted below the steering electronics. A view 
of the octocopter drone with the antenna mount and transmitting antenna is shown on the left in 
Figure ^ The drone had been programmed in such a way that it flew to several pre-defined points 
at a certain distance to the LBA. The position of the drone was measured via a differential GPS, 
which enables a precise measurement of better than 30 cm of easting, northing and height above 
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ground every 0.5 s. The maximum height achieved in the campaign was 52 m above ground. This 
translates into a positioning uncertainty of less than 0.2 ° with respect to the arrival direction of 
the signal. The rotation of the drone with respect to its yaw, pitch and roll axis has been measured 
and corrected for, which allows us to use the same coordinate system for both antennas. A detailed 
description of the octocopter drone including the electronic set-up can be found in Chapter 7 of 

m- 


3.1.2 Reference antenna - BBOC 9217 Biconical antenna 

The BBOC 9217 antenna is an open biconical antenna developed by Schwarzbeck Mess-Elektronik 
0. This custom-built model uses a 4:1 balun which provides a smoother frequency response and 
less environmental dependence. The antenna has a mass of 300 g, a rod length of 25 cm, and 
an opening angle of 53°. As the flight time of the octocopter depends strongly on the overall 
payload, the low mass makes the antenna well suited for the calibration campaign. A photograph 
of the instrument is shown in Figure ^ The antenna has been calibrated by the manufacturer for a 
frequency range between 30 and 1000 MHz. The frequency-dependent isotropic gain of the antenna 
exhibits a resonance frequency of fres = (360 ±5) MHz, which is outside the frequency range of 
the LB A. The systematic measurement uncertainty of the isotropic gain accounts to ±1 dB. The 
reference gain curve can be found in Chapter 5 of 

3.2 Full signal-chain calibration 

The set-up as described above does not allow for an absolute calibration of the LOFAR system. 
The calibration loop has to be connected to the footpoint of the antenna and therefore does not 
include the influence of cables, bandpass filters on the RCUs and the digitization step. Therefore, 
an additional calibration has been performed, using a transmitting antenna that is mounted to a 
crane. The measurement has been performed using the same hardware as QTEl], which allows us 
to also compare results between the different air shower experiments Tunka-Rex [1T3|1 , LOPES 
llTSll and LOFAR. The drone calibration allows for a fast, high and very flexible positioning of 
the transmitting source around the receiving antenna, which provides measurements of the antenna 
directionality. It, however, has the draw-back that the position is not always completely stable. 
The calibration with a crane is more stable, but measurements of the direction dependence are less 
flexible and certain positions may not be reachable. 

3.2.1 Experimental set-up 

The measurements with the reference source were carried out at the LOFAR superteip during a 
two-day campaign in May 2014. As shown in Figures 0 and a crane has been positioned in the 
antenna field of a LOFAR station close to the superterp. A wooden extension of 5.3 meters was 
added to the crane in order to avoid reflections from the metal construction. At the end of this 
extension the reference source and a differential GPS were attached. With this construction the 
reference antenna was positioned at a maximum distance of r = 12.65 ±0.25 m vertically above 
one dedicated antenna. The alignment of the reference antenna with a LBA dipole arm was possible 
with the help of two strings attached to the mount of the reference antenna as shown on the right of 
Figure |8|. 
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Figure 7. Schematic illustration of the experimental set-up using the reference source. The source is sus¬ 
pended from a crane at about 12 m above the chosen antenna. The signal is received with the LOFAR LB A 
antennas and filtered and digitized at the receiver units (RCUs). The data of all antennas of a LOFAR station 
are read out via the LOFAR system using the transient buffer boards (TBB) as it is done for cosmic ray 
measurements. 



Figure 8. Experimental set-up of the measurement with the reference source attached to a crane. Left: 
Crane with the wooden construction and an LBA underneath. Right: Transmitting antenna as mounted on 
the wooden construction in front of several LB As. 


For data-acquisition the LOFAR system has been used. The TBB ring-buffers of the supert- 
erp stations were read out at least 5 times per position. The final data sample, after quality cuts, 
consists of four read-outs, containing 10 ms of data for each of the 48 antennas at a distance of 
r = 12.65 ± 0.25 m vertically above the central antenna. In the measured LOFAR station the same 
configuration is used as during air shower measurements, which ensures that the calibration in¬ 
cludes the full signal chain. 
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Figure 9. Uncalibrated power spectrum as recorded with the LOFAR system for the LBA directly under¬ 
neath the transmitting antenna. The power in the emitted frequencies is at least three orders of magnitude 
larger than the noise background of the system and the Galaxy. Single RFI lines are visible in the band of 
30 — 80 MHz. Their powers are in all cases at least a factor of 100 smaller than the power emitted by the 
calibration antenna and they are located at different frequencies, which makes a confusion of lines unlikely. 


3.2.2 Reference source - VSQ 1000 

The reference source is a commercial product developed by Schaffner, Augsburg in Germany (now 
TESEQ). It is delivered as a combination of a signal generator RSG 1000 and the biconical antenna 
DPA 4000. The RSG 1000 is a comb-generator, generating a spectrum of single frequencies at 
multiples of 1 MHz in the range of 1 MHz to 1 GHz. In the relevant range of 30 — 80 MHz it 
delivers a mean power of 1 /rW per single frequency. It is battery-operated, which makes it ideal for 
measurements in the field. The DPA 4000 biconical antenna is linearly polarized and its directivity 
close to the main lobe follows roughly a cosine squared. This means that misalignments of less 
than 5° with the receiving antenna result in losses of less than 1%. The VSQ 1000 setup is certified 
for fhe 30 fo 1000 MHz frequency range in fhe forward direcfion. An example of a typical VSQ- 
generafed specfrum defecfed by fhe EBA is depicfed in Eigure 


4. Verifying the antenna model 


If a transmitting antenna and a receiving antenna are a loop where the impedances are perfectly 
matched to each other (Eigure the gain of the transmitting antenna can be described as 


Gt{v,d,<p) 


47rr^5'(v,0,^) 

W) 


(4.1) 


where S is the transmitted power density and Pt the transmitted power. The gain G{v,6,<p) is 
proportional to the square of the VEE, \H {v, 6, of the antenna system. Eor a perfectly matched 
system, as it is used in the flying source calibration, both can be translated into each other [CTl - The 
full EOEAR antenna system, however, includes an intentional mismatch to broaden the bandwidth 
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Figure 10. Gain spectra as measured with an LBA in combination with the network analyzer in a loop. 
The spectra are shown for two zenith angle bins for the complete LBA band of 10 —95 MHz. The left 
figure shows 24 spectra for 0 < 10° and the right figures contains 12 spectra in the bin of 40° < 0 < 50°. 
All measurements have been taken at the same azimuth angle. The lines indicate the gain values of the 
antenna model for the corresponding zenith angles of 5° and 45°. Additional filtering outside of the band of 
30 — 80 MHz is not taken into account in the current model and is indicated by the hashed regions. 

The Friis Transmission Equation connects a transmitting antenna, which sends a signal, with 
a receiving antenna detecting this signal [031] in the far-held and is used to determine the gain of an 
antenna. Taking into account the effective area of the antenna, the power at the receiving antenna, 
Pr can be determined by 

Pr(v) = (^) G,(v,e,(^)G,(v,e,0)p,(v) |d,-d,|2, (4.2) 

where Gr and Gt describe the gain of the receiving and transmitting antenna, respectively. The 
last factor \ar-at\^ accounts for any possible mismatch between the polarization of the impinging 
wave and the polarization properties of the receiving device. If both antenna devices match in 
polarization and rehection, it holds that = 1. If the two antennas are aligned, the gain of 

the transmitter Gt is only frequency-dependent. If then Gt of the transmitting antenna is known, 
the Gr can be easily calculated in units of decibel 

Gr,dB{V,0,(p) = 20\ogiQ{47lr) —lOlogiQ^ +Pr,dB - Pt,dB -Gt,dB{v)- (4.3) 

This equation does not include inhuences from ground conditions or rehections on nearby anten¬ 
nas. Based on these equations and the flexible positioning of the octocopter, the directional and 
frequency dependent behavior of the antenna model can be tested. 

4.1 Frequency behavior 

Figure 0 shows a number of frequency dependent gain measurements for two different zenith 
angle bins together with the average expectations given by the antenna model. The strongest signals 
are measured at the resonance of the system near 61 MHz. The figure also shows that the shape of 
the antenna sensitivity remains similar for different zenith bins, while decreasing in gain of about 
5 dBi. 
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Comparing the antenna model in detail with the measurements, it has to be noted that the 
simulated resonance frequency does not fully match the measured position. Foremost, outside of 
the central band of 30 — 80 MHz the quality of the match is significantly worse. This is partly 
expected due to the missing modeling of the filters on the RCU. As those filters are also selectable 
and the frequency region is currently not used for air shower measurements, these mismatches are 
not investigated in detail for this analysis. Inside the central region the differences observed in the 
shape of the distributions are less severe, however, significant. Both issues give an indication that 
the simplified model of fhe anfenna and fhe elecfronics chain fhaf is used for fhe modeling for fhe 
anfenna sysfem mighf nof be fully sufficienf. A renewed simulafion of fhe LOFAR anfennas is fore¬ 
seen in fhe fufure. This will involve a lof more complexify in modeling fhe individual componenfs 
and requires a significanf efforf. For now and given fhe complexify involved, fhe currenf anfenna 
model is kepf and second order correcfions fo fhis model will be derived as parf of fhis analysis. 

Furfhermore, if should be nofed fhaf fhe calibrafion measuremenf is also only a snapshof of one 
day and fime-dependenf changes are hard fo observe. In fhe complefe cosmic-ray dafa-sef, if can be 
observed fhaf fhe resonance frequency shows small shiffs in fhe order of up fo 1 MHz. If has been 
argued fhaf fhe loop in fhe ferminafing wires of fhe LBA, as well as fhe ground plafe accumulafe 
wafer droplefs. This can acf as an additional impedance and change fhe anfenna behavior. If 
has, for example, been observed direcfly fhaf fhe resonance frequency shiffs abouf 0.5 MHz fo 
lower frequencies after parficularly heavy rain. These effecfs will always have fo be considered as 
addifional uncerfainfies. 

4.2 Directional behavior 

The geomefrical sensifivify of fhe anfenna has also been fesfed wifh fhe same sef-up. The frans- 
mission anfenna has been aligned wifh one anfenna arm af differenf zenifh angles, which resulfs 
for fhe emission af 60 MHz in fhe pattern shown in Figure O- The red line is fhe prediction of 
fhe anfenna model wifh average ground conditions. Uncerfainfies on fhe dafa stem mosfly from fhe 
misalignmenf of fhe fransmiffing anfenna and fhe LBA and fhe unsfable positioning. 

The simulafed anfenna model describes fhe direcfional behavior adequafely. This is also frue 
for fhe ofher frequency bins. In order fo esfimafe fhe influence of differenf ground conditions, sim- 
ulafions wifh various values of permiffivify Sr and conducfivify a were performed. These demon- 
sfrafed only a small influence (< 5%) on fhe measured gain per frequency. However, noficeable 
changes (5 — 10%) on fhe shape of fhe direcfional dependence of fhe anfenna pattern were observed. 

Influences due fo humidify or wind are much more difficulf fo esfimafe. Since fhe measure¬ 
menf campaign only delivers a snapshof of fhe condifions as dafa were gafhered af only one day 
wifh a single anfenna, fhe campaign would need fo be repeated fo measure fhe influences, or fhe 
uncerfainfies have fo be estimated ofherwise. 

4.3 Intermediate conclusions 

From fhe measuremenf of fhe anfenna gain pattern, if can be concluded fhaf fhe frequency behavior 
of fhe anfenna model does nof fully mafch fhe observed behavior. The fall-off from fhe resonance 
is differenf in shape in fhe model fhan whaf was measured during fhe campaign. The direcfional 
behavior is on average well described. Reasons could be eifher incomplete modeling of fhe elec¬ 
fronics componenfs or exfernal influences. If should be nofed fhaf varying ground condifions or 
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Figure 11. Gain factor of the antenna as a function of the zenith angle of the incoming signal orthogonal 
to the antenna plane at 60 MHz. The black points show the measurements determined from two calibration 
flights. The red line indicates the predicted gain of the antenna model at 60 MHz. 

humidity do influence the modeled antenna response, especially with respect to the angular depen¬ 
dence. However, no combination of parameters could be found that showed an improved fit to all 
data, which leads us to conclude that incomplete modeling of the system is the more likely cause 
of the discrepancies. 

The calibration using the closed loop is not suitable for an absolute calibration, as it does not 
take the complete system into account. The obtained results, however, give confidence to use the 
current antenna model as basis and to derive additional frequency corrections based on the absolute 
calibration, which is described in the following sections. These corrections will have to account 
for possible frequency dependent damping of the cables and in the amplifiers and bandpass filters 
on the RCUs, as well as imperfections in the antenna model itself. 

5. Absolute Calibration of the system 

For the absolute calibration two different approaches are used: one involving the calibrated refer¬ 
ence antenna attached to a crane and the other a calibration on the diffuse Galactic radio emission, 
which is the dominant noise source. Both approaches are complete end-to-end calibrations, thereby 
including the entire receiver system, as well as cables and electronics. 

5.1 Results of the full signal-chain calibration using a reference source 

From the reference antenna (see section E3) a calibration curve can be obtained that converts the 
voltage traces measured with the LOFAR system to physical units. 

5.1.1 Calibration curve 

The frequency-dependent calibration factor X(v) signifies the translation between the expected 
power Pe{v) in physical units and the measured power Pm{v) in system units. The measured power 
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Figure 12. Calibration factor X for the amplitude as a function of frequency across the LOFAR band as 
derived from the reference source calibration. The dark region denotes the statistical uncertainties of the 
method, while the lighter region illustrates the systematic uncertainties on the absolute scale. 


in each frequency-bin is obtained from the Fourier transform of the measured analog-digital con¬ 
verter (AD) units as The expected power is calculated as the square of the expected 

voltage, divided by the vacuum impedance Zq and combined with the antenna VEL H{v), so that 
for A(v) 


1 |F(v)|2 _ 1 \E{v)-H{v)\^ 


(5.1) 

(5.2) 


In the case of the reference source calibration, E{v,6 = 0,0 =0) is the electromagnetic field 
emitted by the source, the amplitude of which is obtained from the manufacturer H{v) cor¬ 
responds to the VEL from this same direction. Since the source antenna is linearly polarized, only 
the Jxe component of H{v) contributes. Note that A(v)^ is proportional to power, so that A(v) is 
only proportional to the amplitude. 

Eor the analysis, data are used with block sizes of 65400 samples of 5 ns, corresponding 
to a frequency resolution of ~3 kHz in the 1 — 100 MHz range. The bandwidth is restricted to 
30 — 80 MHz. A Gaussian smearing of the edges affecting 5 frequency bins at both ends of the 
spectrum has been applied to the filter to reduce sharp cut-off effects. Signal peaks from the comb 
generator have a width of less than 9 kHz, corresponding to at most 3 frequency bins with this 
resolution. The background noise, as well as single narrowband noise-lines are at least three orders 
of magnitude lower than the signal, and therefore contribute less than 1 % in power. 

The resulting calibration factor A(v) is shown in Eigure |T^ The curve is relatively constant 
over most of the frequency range, with a clear steepening at the upper range of the spectrum. 
This suggests that the antenna response is flatter at the highest frequencies than what is currently 
modeled. This confirms fhe finding of fhe anfenna gain measuremenf, as shown in Eigure [m Ad- 
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Figure 13. Integrated uncalibrated power received in different frequency ranges (markers) as a function of 
distance between reference and receiving antenna. Fitted to the data (full bandwidth) is the power fall-off 
for either a dominating inductive field (dashed line), or for a dominating far held (solid line). 

ditionally, a clear peak is visible around 59 MHz. This is likely caused by the already discussed 
shift in resonance frequency between data and antenna model. The measured peak shifts in corre¬ 
lation with environmental changes, such as rain, and it is therefore rather challenging to match this 
correctly in a model. Different measurements conducted on the same day show stable results with 
variations of about ~1%. Uncertainties are however dominated by systematic uncertainties in the 
emitted power of the VSQ. 

5.1.2 Uncertainties and cross-checks on calibration 

Several uncertainties have to be considered and are given, in the following, as relative uncertainty 
.r on X(v), meaning that the la confidence-interval is described by [X(v)/(I +.r),X(v) x (1 +.r)]. 

The largest contribution to the uncertainties follows from the source antenna, which is dom¬ 
inated by a systematic uncertainty of 16% on the emitted electric field sfrength as given by fhe 
manufaclurer. Signal-insfabilifies caused by flucfuafions in femperafure infroduce anofher sysfem- 
afic 6% on fhe emifted elecfric field sfrengfh QTEl]. This confribufion is expecfed nof fo be presenf 
during fhe measuremenf due fo fhe same fype of weafher fhroughouf fhe campaign, buf will affecf 
fhe applicabilify of fhe calibration af ofher environmenfal condifions. Addifionally, possible mis¬ 
alignments between the source and the receiving antenna influence the result. Offsets are estimated 
to be less than 0.25 m in altitude and 5° in rotation, resulting in uncertainties on the calibration 
factor ^(v) of up to 5%. Considering all influences of the measurement setup, such as antenna 
cross-talk and signal reflections on nearby surroundings, is believed to yield uncertainties of no 
more than 5%. 

For an accurate approximation of the signal, the source needs to be located in the far field 
where fhe power P falls off as P oc as discussed in secfion ||. The determining distance for our 
case is r >> A, where A = 10 m for 30 MHz. This questions whether our measurements are truly in 
the far-field. However, the integrated measured powers reveal a good fit to a far field approximation 
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Uncertainty (a) 



Value [%] 

Antenna-by-antenna 


Variations between antennas 

1 



Total 

1 

Event-by-event 


Environmental (excl. source) 

5 



Total 

5 

Calibration 

Method-specific 

Source emission 

16 



Far-field 

5 



Total 

17 



Set-up 

5 


Campaign-specific 

Alignment 

5 



Source temperature stability 

6 



Measurement variations 

1 



Total 

9 


Table 1. Summary of the uncertainties on the calibration curve in amplitude that have to be considered for 
the full signal-chain calibration using the reference antenna. 


for all frequencies even at distances r < A as shown in Figure |TJ. While this suggests the correct 
use of a far field approach, spill-over effects from the near field cannot be excluded completely. 
Therefore, an additional uncertainty originating from a far-held approximation is approximated to 
be 0.5 dB (~ 6% in amplitude) JTTl]. 

During the calibration measurement the additional 47 dual-polarized antennas were also read 
out in the LOFAR station. This could allow one to cross-check the calibration value and check 
for uncertainties between antennas. However, in order to be able to do so, one needs to make 
an assumption about the emission pattern of the calibration source since the manufacturer only 
provides values for the main lobe. This additional antenna model increases the complexity of the 
test. Differences between the expectation and the measurement can either be attributed to issues 
with the model of the LBA or the one of the source antenna and neither can be excluded. As a hrst 
test indeed shows an unexpected but smooth trend in the calibration curves, we have refrained from 
using receiving dipoles not directly under the reference antenna. 

When comparing the integrated power between different measurements on the same day, the 
statistical fluctuations are less than 1%. From the Galactic calibration (see section which can 
take into account all antennas, we have estimated that the variation between antennas is about 1%, 
if the antenna is fully functioning, which is the case for the measured antenna. Additionally, also 
following the Galactic calibration, environmental effects on the LBAs are determined to contribute 
5% uncertainties. An overview of all different contributing uncertainties can be found in Table |T]. 
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Using conservative estimates the total uncertainties sum up to 


^X(v) 

—= ± 1% (antenna-by-antenna) ±5% (event-by-event) (5.3) 

X(v) 

± 17% (method-specific) ±9% (campaign-specific). (5.4) 

More details about this analysis can be found in [051] . 

5.2 Results of the calibration using the Galactic emission 

The LBAs have a large field of view and data collected with the system are known to be sky noise 
rather than system noise dominated. With the cosmic-ray data collected since 2011 a precise and 
continuous reference of the noise floor is available which is used for calibration purposes. 

The measured noise is a combination of sky and system noise. The system noise temperature 
Tsys{v) is defined as the sum of the receiver noise temperature Trec{v) and the electronic noise 
temperature Tgiedv): 

Lysiv) = Treciv) +Telec{v). (5.5) 

The latter combines all noise contributions from filters, amplifiers, cables and receivers. In the 
absence of RFI, the receiver noise temperature only contains emission from the Galaxy, so that 
Trec{v) = Tsky{v) The combined noise floor of electronics and sky, Tsys, is compared to 

predictions from models of the diffuse emission of the Galaxy in order to obtain an absolute cali¬ 
bration. 

5.2.1 Predictions of the Galactic emission 

In this analysis, two different models are used for the prediction of Galactic radio emission; LFmap 
[0711 and the Galactic Sky Model (GSM) [^Bl] . Both work by interpolating published reference sky- 
maps in the tens of MHz to several GHz range to any desired frequency in that range, and calculate 
a sky map Tsky{v,a,8) with (a, 5) being equatorial coordinates. The two models differ by less 
than 5% in the LB A bandwidth and quote similar uncertainties. Two types of uncertainties have to 
be distinguished. Both models fit the reference maps, which means that there is a contribution of 
the quality of the fit and of the absolute scale that the fits inherit from the original measurements. 
The fit is reported to describe the underlying maps with an accuracy of usually 5% and only for 
some maps with up to 10% across the whole fitted range of frequencies. 

The two models directly inherit the absolute scale of the reference maps and can be only as 
good as the data. The highest quoted uncertainty of the original data (20% in power) [^El , ^ is 
therefore assumed as a conservative estimate for the overall scale uncertainty. This conservative 
estimate also accounts for the fact that the LOFAR bandwidth has been less widely covered in 
surveys and some reference maps are incomplete for a few regions with high emission. Future 
dedicated work might be able to reduce these uncertainties to less conservative values. 

As differences between LFmap and GSM are relatively small, only the LFmap model of the 
Galactic radiation is used for the analysis. An example of a map of the emission at 60 MHz as 
obtained from LFmap is shown in Figure [141 . 
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Figure 14. Map of the Galactic radio emission, generated by LFmap [P7[] at 60 MHz. The upper region is 
visible from the location of LOFAR at 17.3 hours Local Sidereal Time (LST). The black line corresponds to 
the horizon and the dashed region is not visible at this specific time. 


The power expected from Galactic radio emission can be calculated via the Rayleigh-Jeans 
law. The received spectral power Pv,z in one of the LOFAR dipoles (X) is given by 


2kB 


Pv.x = ^v^ T{v,e,(i))\Hiv,e,^)\Ua 

Ja 

<Pv,x> = [ T{v,e,(l)){\Jxe\^ + \Jx^\^)da. 

Ja 


(5.6) 

(5.7) 


using equation More precisely, equation follows as a result of a time-averaged antenna 
response, determined by the received voltage Vx in a dipole X as the response to unpolarized waves 

<V'i> = <{E{v,e,^)-Hx{v,d,^)f> (5.8) 

= < {JxQ sm{(o)f + {Jx(j, cos{(o)f + 2 -Jx9JxiP sin(a)) cos(a)) > (5.9) 

= 2^\Jx9\'^ + \Jx(j>\^)- (5.10) 

The antenna-based coordinates (0,0) are calculated from (cc,d) using the Local Sidereal Time, 
and the location of the LOFAR superterp at latitude 52.92 °N and longitude 6.87 °E. As the visible 
sky region changes over time, the power received from the Galaxy is also time-dependent. If such 
a variation is observed in the data, the system is sensitive to and likely dominated by the sky noise. 

5.2.2 Contribution of the electronic noise 

The measured power is depicted in Figure [T^ as a function of Local Sidereal Time. A clear vari¬ 
ation is visible. In order to match simulated powers from the Galactic emission to the detected 
data, the contribution from the electronic noise needs to be included. The system noise has been 
measured during deployment and commissioning PD| , ES 0- those measurements, however, 
only intended to show that the system is sky noise dominated and did not focus on providing an 
uncertainty. Thus, our data are used to directly determine the most probable fraction of sky noise. 


- 19 - 






Figure 15. Integrated median uncalibrated noise power, as a function of the Local Sidereal Time, for the 
50 — 60 MHz sub-band. Also shown is the predicted received power in both dipoles (dashed red lines), 
before (left) and after (right) applying electronic noise corrections. 



Frequency [MFIz] 


Figure 16. Calculated values for the contribution of the sky temperature to the total system temperature. 
Shown are both the results from this analysis (blue lines with shaded uncertainty) and an approximation 
of the median values derived from the observation of 3C295 [0 (dashed dotted red line). The dashed line 
indicates a fraction of 0.5. 


The most probable noise offset due to the electronic noise has been found using a least-squares 
fit. The data was binned in LST-bins of about 15 minutes. With the current complete data-set 
every bin then contains the noise background of an average 34 air showers per antenna. As the 
electronic noise is frequency-dependent, offsets need to be determined in frequency sub-bands. 
Here, sub-bands of 10 MHz (i.e. ranging [30 — 39,40 — 49,... ,70 — 79] MHz) were chosen. The 
reduced was calculated for every combination of electronic noise offset per frequency band 
and the simulations of the Galaxy with respect to the binned data. The uncertainties per sub-band 
correspond to noise corrections at the point where x^ = Xmin + 1’ which is a conservative estimate 
of the uncertainty. 
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With these values, a measure for Tsky/Tsys can be constructed per sub-band Vo as 

^Vq, sky _ _ ^Vq, sky _ ^ 

Tvn,SyS ( ^VQ: noise , rp 

\ ti j“rrvo,iiy 

*'30 — 80. mean 

Here, the temperature sky in a certain sub-band Vq is determined as the average sub-band voltage 
divided by the average VEL amplitude of the antenna model in each sub-band. The electtonic 
noise offset per band is divided by the average value of the VEL in the full 30 — 80 MHz range 
for normalization. Resulting values are depicted in Eigure |T^. Values are comparable to what has 
been established earlier [0. The largest discrepancy is near the resonance frequency, which is most 
easily affected by using a slightly different antenna model. As the uncertainties on the astronomical 
method can no longer be obtained, no significance of this discrepancy can be given. 

The variation of both measured spectral power and predicted power before and after noise 
corrections are shown in Eigure |T^ for the 50 — 60 MHz sub-band, where LOEAR is most sensitive 
due to the resonance frequency of the dipole. The figure shows that the predicted curves match 
the data better after the noise correction. Subtle systematic differences indicate that an improved 
understanding of the contribution of the electronic noise and the antenna model can positively affect 
this analysis. A similar behavior can be observed for the other frequency bands. The obtained 
noise corrections have been added to the simulations of the sky noise for all further analyses. It 
is foreseen to measure the noise contribution in a dedicated measurement campaign to reduce the 
uncertainties. 

5.2.3 Calibration curve from comparison to data 

In a way identical to the reference source calibration (eq. a calibration curve for the full signal- 
chain can be established for the predicted Galactic emission using the same antenna model. The 
expected power is composed of both Galactic radiation and electronic noise combined, while the 
measured power in each frequency bin is again calculated as |.^(v)p. 

The resulting calibration curve is shown in Eigure [T^ This method predicts an increase in 
calibration factor towards higher frequencies. Noticeable bumps are visible at 50, 60 and 70 MHz, 
mostly as a result of binned determination of the electronic noise contribution. This indicates that a 
different measurement of the electronic noise can be used to improve this analysis. This, however, 
requires a potentially time-consuming dedicated measurement effort. Alternatively, the current 
results could be smoothed before applying the calibration. 

5.2.4 Uncertainties 

The largest source of uncertainties on A(v) follows from the contribution of the electronic noise to 
Tsky/Tsys- This uncertainty propagates into a systematic uncertainty of 37% on the calibration curve. 
It should be noted that this is a conservative estimate of the uncertainty. The second important 
contribution is caused by the prediction of the Galactic radiation. The absolute scaling uncertainty 
of different reference maps ranges up to 20% in temperature, while the use of different models 
yields a difference of 5%. The scaling uncertainties in temperature result in uncertainties on A(v) 
of 9% and 2%, respectively. 
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Figure 17. Calibration factor for the measured amplitude as a function of the frequency as derived from 
diffuse Galactic emission. The dark region denotes the statistical uncertainties of the method, while the 
lighter region illustrates the systematic uncertainties on the absolute scale. 


The statistical uncertainties contain several aspects. Per data sample, they contain all used 
dipoles (up to a maximum of 96), as well as data collected over the course of several years. There¬ 
fore, all occurring weather circumstances (ranging from sub-zero up to 35 °C, from dry conditions 
to rain or snow) are included in the spread of the calibration curve. Measured statistical fluctua¬ 
tions have a standard deviation of up to 5% over the full frequency range. Fluctuations are of the 
same order when analyzing all events for any single antenna, so that it can be concluded that this is 
mostly dominated by environmental circumstances rather than fluctuations between antennas. An 
overview of all different contributing uncertainties can be found in Table The total uncertainties 
sum up to 

= ±1% (antenna-by-antenna) ± 5% (event-by-event) ± 38% (calibration). (5.12) 


5.3 Comparison of calibration methods 

The two methods to obtain an absolute calibration are compared to each other and the consistency 
with other experiments as reported in the literature is checked. Since the absolute signal amplitude 
scales with the energy of the air shower (e.g. most recently the consistency between ex¬ 

periments is especially relevant for studies of cosmic ray properties and the corresponding energy 
scale. 

5.3.1 Comparison of Galactic and reference source calibration 

Both calibration methods result in frequency-dependent calibration factors with a very similar be¬ 
havior as a function of frequency (Figure |^. The two curves are fully compatible below 60 MHz, 
and only above this frequency slight deviations are visible. Here, it is interesting to note that the 
shape of the two curves also deviates. This can probably be attributed to the fact that the two 
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Uncertainty (a) 


Value [%] 

Antenna-by-antenna 

Variations between antennas 

1 


Total 

1 

Event-by-event 

Environmental 

5 


Total 

5 

Calibration 

Choice of sky model 

2 


Absolute scaling of model 

9 


Relative scaling of model 

5 


Electronic noise 

37 


Total 

38 


Table 2. Summary of the uncertainties on the calibration curve in amplitude that have to be considered for 
the calibration on the diffuse emission from the Galaxy. 



Figure 18. Calibration factors as function of frequency across the LOFAR band for Galactic and reference 
source calibration. Both calibration curves contain statistical uncertainties of the method in the dark region, 
with systematic uncertainties illustrated by the lighter region (dashed for Galactic, hlled for terrestrial). 


methods use different types of signals. While the reference source calibration exploits signals of 
several orders of magnitude above the noise level, the Galactic calibration relies on the noise level 
itself. Being essentially a simple dipole, the LBAs are mostly sensitive to the resonance frequency, 
meaning that for higher frequencies the antenna becomes too long (inductive) and its impedance is 
no longer small with respect to the LNA. Thus, the gain of the LNA decreases and the contributions 
of the noise budget accumulated in the coaxial cables and the several amplification stages becomes 
relevant. This, however, does not affect the strong signals of several orders above the noise level. 
Consequently, the two curves show a slightly different shape with respect to the antenna model that 
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does not yet include the above mentioned cables and amplifiers. 

Furthermore, the Galaxy is not a point-source like the reference antenna. The measurement 
is always an integration over all visible angles. This means that all discrepancies between antenna 
model and the antenna behavior with fold into the measurement. This is different for the reference 
source, where only the directions directly above the antenna are calibrated. The two components 
of the antenna model, directivity and frequency dependence, can therefore not be decoupled with 
the calibration on the Galactic emission. 

Both curves show a steepening slope towards the higher frequencies, which confirms the an¬ 
tenna calibration measurements with the drone that showed that the current antenna model de¬ 
scribes a steeper fall-off than seen in data. Potential reasons are manifold. For example, coax 
cables are known to yield a signal-attenuation that increases with frequency, which has not been 
considered in the current model. Simplified corrections improve the situation, but this cannot fully 
explain the observed steepening. The obtained electronic noise corrections for the Galactic cali¬ 
bration result in visible bumps around 50, 60 and 70 MHz. Finally, signal reflections and known 
cross-talk of closely spaced antennas can affect the averaged Galactic measurements. 

The reference source calibration is essentially only a snapshot of the data, which means that 
time-varying influences had to be estimated rather than measured. For the Galactic calibration, a 
better time-average is possible, however, the electronic noise introduces unfortunate features and 
currently a rather large scale uncertainty. Both results are currently completely compatible and of 
similar quality. 

5.3.2 Comparison to other experiments and previous calibrations 

A calibration using the same source antenna has been performed on radio antennas at both LOPES 
and Tunka-Rex. The measured voltage traces after the application of the appropriate calibration 
curves at the three different experiments is shown in Figure [T^ As the calibration technique is 
designed to calibrate the detected pulses to values reported by the manufacturer of the VSQ 1000, 
the pulses from different experiments are required to be of the same amplitude after normalization 
to the same source-distance (10 m). Differences in a single measurement only occur due to different 
hardware response, such as group delays and sampling frequencies. The figure shows that the 
LOPES pulse is slightly lower than the ones from the other experiments. It should be noted that 
the LOPES pulse is an average of a large number of measurement campaigns and thereby slightly 
different. Since it is compatible within the campaign uncertainties, it gives an indication that the 
remaining systematic uncertainties between the experiments are well understood. 

It should be noted that the reference antenna has been re-calibrated between the published 
results in QTEl ] and this campaign. The newest calibration has been performed in accordance with 
free-space circumstances, which matches the needs of air shower experiments. Differences be¬ 
tween the old reference values (reflective ground) and the new values are typically a factor of two 
in amplitude [^31] . Eurthermore, the magnitude of the absolute uncertainty has been improved 
considerably. 

Having the same reference source at the three different experiments, allows us to ignore some 
absolute scaling uncertainty when comparing results from the different experiments in future pub¬ 
lications. The largest contribution to the scale uncertainty as shown in Table |T] are the method- 
specific factors. As opposed to the campaign specific ones, these uncertainties fully cancel out. 
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Figure 19. Electric field as function of time for the three calibration measurements at (from left to right) 
Tunka-Rex, LOFAR and LOPES, with amplitudes scaled to match a source at 10 meter distance, up-sampled 
and hltered to a common band of 43 — 74 MHz. The dashed lines represent the calculated Hilbert envelopes 
m for the three pulses. They allow for a better comparison of the data that are all taken with different 
sampling frequencies. The size of the systematic campaign specific uncertainty on the maximum amplitude 
is indicated with the grey bar. These uncertainties are expected to be of the same order of magnitude for 
Tunka-Rex and LOPES. The pulses are offset in time from each other for a better visibility. 


when comparing air shower measurements at LOPES, LOFAR and Tunka-Rex. This means that 
the reconstructed energy of cosmic rays from radio emission of air showers becomes directly com¬ 
parable, with small systematic uncertainties. 

6. Cross-check with CoREAS air shower simulations 

The calibration values as obtained above have been compared to predictions from air shower simu¬ 
lations. The air shower simulations predict the pulse strength delivered per ideal antenna, given that 
shower axis position, energy, height of shower maximum and arrival direction are known. When 
folding the ideal signal strength with the specific antenna model, a direct comparison of model and 
data is possible m- 

A selection of high-quality air showers is routinely simulated with CoREAS |@] to determine 
the height of the shower maximum. In order to match the absolute signal strength to the (still) 
uncalibrated LOFAR data, a scaling factor S is introduced, which compensates for the missing 
calibration. 

We compare the simulated pulse energy density (CoREAS and antenna model) and the mea¬ 
sured pulse energy density. Each energy density is obtained by integrating the squared pulses over 
/ = 11 samples of 5 ns each around the pulse maximum. This integrated value is corrected for the 
average noise level in each time-bin by subtracting its contribution from the integrated value [HDH . 
To obtain the total energy density (/,„), both instrumental polarizations are added. The scaling fac¬ 
tor S is the ratio of the (uncalibrated) energy density derived from measurements (/„,) and the same 
quantity obtained from CoREAS simulations in combination with the antenna model, (4): 
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Figure 20. Scaling factor S as obtained in the time domain for the air shower simulations (CoREAS) 
and the two calibration methods (Galaxy and Reference source). Depicted values and uncertainties are 
the average and its standard deviation obtained from different events. The colored bands correspond to 
systematic uncertainties on the used method as discussed in the corresponding sections. 


S=Y [ 5 ]=£^. ( 6 . 1 ) 

In future work, we plan to implement the two calibration curves in the data analysis pipeline 
and to perform an in-depth study of the agreement of air shower simulations with data, including 
the timing, polarization and frequency behavior. Here, as a first step, we compare the scaling 
factors S obtained from every method. 

As shown in 0, the energy used as input to the air shower simulations has to be correct to 
obtain a stable scaling factor. As the signal amplitude scales directly with the energy of the shower, 
the accuracy of the scaling factor S between LOFAR data and CoREAS simulations is severely 
determined by the energy resolution and energy scale of the particle detectors. Furthermore, the 
method used to combine the fitting of particle and radio data determines the scaling factor. 

The factor S describing the same conversion of integrated powers has been obtained for both 
the reference source calibration and the Galactic calibration. In the case of the reference source, 
measured pulses^ are squared, integrated in the time-domain and compared to the expectations. It 
was found that subtracting the noise-level contributes less than 1% in cases of rather strong signals 
like the ones from the reference source. For the Galactic calibration, emission is not restricted 
to pulses, and the entire trace is integrated and normalized to the same time-width. As power 
is conserved in Fourier transforms, the powers for the Galactic contribution are integrated in the 
frequency-domain. 

We would like to stress that the method by which S is obtained varies significantly from the 
calibration method used earlier in this analysis to obtain X (v). An important advantage of X (v) is 
that it allows for a discussion of frequency-dependent characteristics, which the scalar value S does 

frequency comb corresponds to pulses in the time-domain 
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not. X{v) is therefore the desired quantity. Air shower pulses, however, are usually shorter than 
60 ns, which does not allow us to determine A(v) with a high frequency-resolution. Furthermore, 
it is expected that air shower pulses show a dependency of frequency as function of distance to the 
shower axis and the height of the shower maximum. These effect have yet to be experimentally 
confirmed. Consequently, we cannot derive a third independent A(v), but have chosen to use S, 
which is independent of frequency. 

The resulting values for S for all three methods are depicted in Figure including their statis¬ 
tical and systematic uncertainties. The statistical uncertainties of the scaling value S for CoREAS 
are determined by the energy resolution of the air shower measurements. Systematic uncertain¬ 
ties include the absolute scale of the energy of the measured cosmic rays, and uncertainties due 
to the air shower simulations of CoREAS. Eor the method currently used, the uncertainty on the 
absolute energy scale is estimated conservatively at ~50% [J71],[P^. The precise uncertainties due 
to emission models or hadronic interaction models are more difficult to establish. They can be 
estimated by the comparison of the two air shower simulations CoREAS and ZHAires 0 and how 
they describe the EOEAR data. The results differ by less than 30% in S for a small test-sample of 
showers [^31] . It cannot be established whether the differences are due to the underlying particle air 
shower simulations or the implementation of the calculation of the radio emission. The difference, 
however, provides an estimate of the uncertainty. 

Both, Galactic and reference source calibration result in very similar values for S, which was 
expected due to the agreement shown in Eigure [T^ We have to note that there seems to be some 
tension between CoREAS and the calibration values. However, the observed differences are of 
the order of magnitude of the systematic uncertainties, so the findings are still in agreemenf wifh 
a correcf predicfion of fhe air shower emission sfrengfh wifh CoREAS. A more defailed compari¬ 
son is needed fo draw firm conclusions abouf fhe compafibilify. A comparison also including fhe 
polaiizafion of fhe signal is currenfly being sfudied, buf is beyond fhe scope of fhis paper. 

7. Conclusion and Outlook 

We have presenfed fhree differenl mefhods fhaf have been used fo calibrafe fhe low-band anfennas 
of fhe EOEAR radio felescope. The firsl mefhod provides defailed informalion abouf fhe anfenna 
response. By posifioning a reference source af differenf posifions around an anfenna, ifs direc- 
fionalify as well as fhe frequency response have been be measured and compared fo models of fhe 
hardware. The anfenna model for fhe EOEAR EBAs has shown fo be in reasonable agreemenf wifh 
fhe measuremenfs. If provides fhe foundation for secondary corrections fhaf have been derived 
during fhe absolufe calibrafion performed in fhis analysis. 

Based on fhis anfenna model, fhe fwo ofher presenfed calibrafion mefhods have been used 
for an absolufe calibrafion. Bofh mefhods, a calibrafion on a reference source and one on fhe 
diffuse Galacfic emission, deliver resulfs fhaf are in good agreemenf wifh each ofher. The reference 
source mefhod is subjecf fo a sysfemafic scale uncerfainfy of 19%, while fhe Galacfic calibrafion is 
dominafed by fhe uncerfainfy of fhe noise confribufion of fhe sysfem of 38% in amplifude. 

Bofh measuremenf campaigns demonsfrafed fhaf a calibrafion wifh reference sources requires 
considerable experimenfal efforf and has fo be repeafed several times fo reduce sysfemafic depen¬ 
dencies on weafher and femperafure condifions. Due fo fhe calibrafion of fhe reference sources, fhe 
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largest systematic uncertainties can hardly be reduced. Reference sources allow for very flexible 
testing of different components of the system and provide valuable information about the model 
of the antenna. They are an essential tool to understand whether the response of the hardware has 
been modeled correctly in software. Better calibration measurements in anechoic chambers are 
only possible in very large chambers at these low frequencies and will not include ground effects. 
It is therefore likely that field measurements will remain an important tool. 

Calibrations on astronomical sources or the diffuse emission of the Galaxy are less time- 
consuming. If the experiment is located in an environment with little external noise and if the 
antenna is designed to be sky noise limited, the data to perform this calibration can be gathered 
during regular operations for every single antenna of the experiment. The largest systematic un¬ 
certainty at LOFAR is the missing detailed knowledge of the contribution of the electronic noise. 
Dedicated measurements can reduce these uncertainties within a manageable amount of effort. This 
will bring the systematic uncertainty down to the uncertainty of the sky models, which are better 
than 10% in amplitude and make the Galactic calibration compatible in accuracy to the reference 
source. 

Further improving the antenna model itself, will also help to reduce the uncertainties of the 
calibration curves. Using the campaign data and a new antenna model, the calibration curve can be 
updated without additional experimental effort. Every calibration curve is only valid in combina¬ 
tion with the antenna model. It seems worthwhile to consider to also measure the phase response of 
the LOFAR system, which is included in the antenna model and especially relevant for broad-band 
transient data. 

The calibration curve obtained from the reference source will be used in future studies as an 
absolute reference scale to further understand the radio emission of air showers. Using this scale 
makes the results of LOFAR directly comparable to those from other experiments calibrated with 
the same source since the largest systematic uncertainties will not be relevant in the comparison. 

Despite a desirable improvement, the current limiting uncertainty of 19% is already competi¬ 
tive with respect to air shower physics. The amplitude of air shower pulses scales with the energy 
of the shower. If dependent only on the radio signal, the amplitude calibration will therefore de¬ 
termine the energy scale uncertainty. The best current energy scale uncertainty is delivered by 
Fluorescence telescopes and is quoted at 14% HIS- Particle arrays have to be calibrated against 
other methods, for example air shower simulations, and are therefore usually subject to system¬ 
atic uncertainties larger than 19%. Their decreased sensitivity to the chemical composition of the 
cosmic rays adds another systematic uncertainty. 

At LOFAR, the systematic uncertainty from the energy scale of the particle detector array 
limits conclusions about the correctness of the models of the radio emission of air showers. From 
our analysis, one could argue that there is some tension between the calibration measurements and 
the values derived from the air shower emission as predicted by CoREAS. However, systematic 
uncertainties are substantial, and no significant disagreement was found. A more dedicated com¬ 
parison of air shower models and measurements, also taking the signal timing, frequency response 
and polarization into account, will shed more light on this issue and is currently underway. 

The calibration curve for air shower measurements can easily be applied to other measure¬ 
ments with the Transient-Buffer Boards at LOLAR to obtain, for example, an absolute flux scale 
for pulsars. Collecting suitable data to test this effort is currently ongoing. 
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